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ABSTRACT 



We present the first astronomical detection of a diatomic negative ion, the cyanide anion CN", as well as quantum mechanical 
calculations of the excitation of this anion through collisions with para-H2. CN" is identified through the observation of the J = 
2-1 and J = 3-2 rotational transitions in the C-star envelope IRC +10216 with the IRAM 30-m telescope. The U-shaped line profiles 
indicate that CN", like the large anion CjH", is formed in the outer regions of the envelope. Chemical and excitation model calculations 
suggest that this species forms from the reaction of large carbon anions with N atoms, rather than from the radiative attachment of 
an electron to CN, as is the case for large molecular anions. The unexpectedly large abundance derived for CN", 0.25 % relative to 
CN, makes likely its detection in other astronomical sources. A parallel search for the small anion C2H" remains so far unconclusive, 
despite the previous tentative identification of the / = 1-0 rotational transition. The abundance of C2H" in IRC +10216 is found to be 
vanishingly small, < 0.0014 % relative to C2H. 
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1. Introduction 

The molecular anions detected so far in the interstellar and 
circumstellar gas are all fairly heavy linear carbon chains 
with three or more carbon atoms, and with neutral counter- 
parts with large electron affinities: C4H , CeH , C gH , C^N_, 
and C-sN" (McCarthv et al., 2006; Cernicharo etall 12007 



CsN' ( McCarthy et al., 2006; Cernicharo et all 
'Brunken et al., 2007a; Remiian et al.L l2007t iThaddeus et al 



|2008; .Cernicharo et al.. .2008.) . The abundance of these anions 
relative to the neutral counterparts increases with size and with 
the electron affinity of the neutral molecule, as expected for 
formation by radiative electron attachment (Herbst & Osamur^, 
[2OO8 ). On inspection, however, that process fails to explain the 
abundance of the shortest observed anions, in particular C4H" 
and C3N . In IRC +10216, a carbon star envelope where both 
anions are found, C3N" has an anion-to-neutral abundance ra- 
tio about 50 times higher than that of C4H , indicating that 
other formation processes may be at work {C ernicharo et al., 
2007| ;lThaddeus et al.Ll2008tlAgilndezLl2009l;lc'oi-diner & MiUaii 



2009|). Studying the astronomical abundance of even shorter an- 
ions, in particular C2H" and CN , whose formation by radiative 



* Based on observations carried out with the IRAM 30-meter tele- 
scope. IRAM is supported by INSU/CNRS (France), MPG (Germany) 
and IGN (Spain). 



electron attachment is very slow, should help answer this ques- 
tion. 

In this Letter we describe the identification in IRC +10216 
of CN" and the results on a parallel search for C2H . We also 
present quantum mechanical calculations of the collisional exci- 
tation of CN" by para-H2, using the calculated rate coefficients 
to model the observed lines. The chemistry of CN" in space is 
also briefly discussed. 

2. Observations and identification of CN 

The C2H" and CN" anions are closed-shell molecules whose 
rotat ional spectrum has been recently measured in the labora - 
tory (iBriinken et al.Ll2007bHGo"ttlieb et al.Ll2007HAmanoll2008h . 
Their electric dipole moments are 3.1 and 0.65 Deby e respec- 
tively dBrunken et al.L l2007bl; iBotschwina et al.L Il995h . 

The present astronomical observations were carried out to- 
wards IRC +10216 with the IRAM 30-m telescope on Pico 
Veleta (Spain). The J = 1-0 rotational transition of CN" at 1 12.3 
GHz was observed before 2009 in the course of a /I 3 mm spec- 
tral survey (Cernicharo et al. in preparation). In IRC +10216's 
spectrum, this line is severely blended with a strong component 
of the ^03/2 J = 81/2-79/2 transition of CeH (see Fig. [B- The J 
= 2-1 and J - 3-2 rotational transitions of CN", at 224.5 and 
336.8 GHz respectively, were observed between January and 
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Table 1. Observed line parameters of CN 



Ta (mK) 



Transition 



vo" 
(MHz) 



(MHz) (kms-') (Kkms-') 



J= 1-0 

J = 2-i 
J = 3-2 



112264.8 
224525.1 
336776.4 



112264.8'' 

224525.4(5) 

336777.0(12) 



14.5'' 
14.5'' 
15.0(10) 



- 0.07(3)" 

0.23(7)" 

0.13(2) 



Note.- Number in parentheses are Icr uncertainties in units of the last 
digits. " Frequencies derived from the rotational constants reported by 
lAmanol ((2008). * Vexp is the half width at zero level. ^ Fixed value. " 
Highly uncertain estimate. Line severely blended with a strong CeH 
line. "^ Line blended with a SiC2 V3=2 line. 



April 2010 with the new dual polarization EMIR receivers oper- 
ating in single side band mode. The rejection of the image side 
band was 13-15 dB at 224 GHz and 20-30 dB at 336 GHz, de- 
pending on the polarization, as measured with strong lines. The 
local oscillator was shifted in frequency to identify possible con- 
tamination from the image side band. The backends were two 
autocorrelators with 2 MHz and 320 kHz channel spacings, re- 
spectively. The pointing and focus of the telescope were checked 
every 1-2 hours on Mars and on the nearby quasar OJ 287. To 
obtain flat baselines, the secondary mirror was wobbled by 180" 
at a rate of 0.5 Hz. The zenith sky opacity at 225 GHz was typ- 
ically ~ 0.1, resulting in system temperatures of 140 K at 224 
GHz and of 800 K at 336 GHz. The total integration time per 
polarization was 3 h at 224 GHz and 9.5 h at 336 GHz, yielding 
a rms noise of T^ ~ 2 mK per 2 MHz channel at both frequen- 
cies, after averaging both polarizations. 

Following our initial detection in IRC +10216 of a T* ~ 
3 mK line at the freque ncy of the 7=1-0 transition of C2H" 
(ICernicharo et al.L 12008 *). we searched from January to April 
2010 for the 7 = 2-1 transition at 166.5 GHz. The line was not 
detected with a T^ rms noise level of 0.6 mK per 2 MHz chan- 
nel, casting doubt on the tentative identification of the C2H" 7 = 
1-0 line. 

The CN observed lines are shown in Fig.[T]and the derived 
line parameters are given in Table [T] The J - 3-2 transition of 
CN" is shown in the top panel of Fig. [T] It appears as a U-shaped 
line with the expected half width (Vexp = 15 ± 1 km s ') that 
agrees in frequency to within 0.6 MHz with that of the CN" 
transition. The 7 = 2-1 transition of CN", shown in the middle 
panel of Fig. [T]with a spectral resolution of 2 MHz and of 320 
kHz (2.7 and 0.4 km s"' respectively), coincides with a broad 
spectral feature with a complex shape that is unusual for IRC 
H- 102 16, since it is neither U-shaped, flat-topped or parabolic. 
It is best explained as a blend, as shown in Fig. [T] that can be 
well fitted with two components, one U-shaped with a half width 
Vexp of 14.5 km s"' centered at the frequency of the 7 = 2-1 
transition of CN" (see Table [T]|, the other with a parabolic pro- 



file, a half width Vexp of 15 ± 3 km s" , and a rest frequency of 
224518.3 + 1.5 MHz that is close to that of the 102,9-92,8 rota- 
tional transition of SiC? in the V3=2 vibrational state (224519.7 
MHz: llzuha et all [19941) . Since other V3=2 fines of SiC2 with 
similar intrinsic strengths have similar shapes, half widths (Vgxp 
= 8-15 km s"'), and intensities {T^ ~ 20 mK) in our A 0.9 mm 
data (Kahane et al. in preparation) as our fitted parabolic com- 
ponent, there is little doubt that this component comes from 
SiC2. We note that the CN" 7 = 2-1 transition has several hy- 
perfine components due to the nitrogen quadrupole, which can 
be grouped into three blocks lying at 224523.9, 224525.1, and 
224527.2 M Hz, with relative line strengths of 0.27, 1, and 0.12, 
respectively (JGottlieb et al.L l2007h . Due to the severe blending 
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Fig. 1. Spectra of IRC +10216 covering the 7 = 1-0 to 7 = 3-2 
transitions of CN". Grey horizontal boxes mark their expected 
positions based on the laboratory frequencies and a linewidth of 
29 km s"'. Shaded areas show the fits to the line profiles ob- 
tained with the CLASS method shell. The high spectral resolu- 
tion spectrum of the 7 = 2-1 line shows the expected position 
of the different hyperfine components with their relative intrin- 
sic strengths. The intensity scale is expressed as T^, antenna 
temperature coiTected for atmospheric absorption and antenna 
ohmic and spillover losses. To transform into main beam bright- 
ness temperature (Tmb) in this Figure and in Table [T] divide by 
0.78, 0.65, and 0.44 at 1 12, 224, and 336 GHz, respectively. 



with the SiC2 V3=2 and to the limited sensitivity of the astro- 
nomical observations, only the strongest hyperfine component 
is clearly visible in the spectrum of IRC +10216, while the mid- 
dle strength component is hidden between the two stronger fitted 
lines (see Fig.[T]), and the weakest hyperfine component lies be- 
low the noise level of the spectrum. Finally, the bottom panel of 
Fig. [U shows the spectrum covering the CN" 7=1-0 transition, 
which is heavily blended with a strong line of CeH. The limited 
spectral resolution (1 MHz) and the broadening of this CN" line 
by the hyper fine structure (there a re three components separated 
by 1-2 MHz: lGottfieb et alJl2007l) makes it difficult to determine 
the relative contributions of CgH and CN" to the observed line. 
There are no good candidates other than CN" for the carrier 
of the 336777.0 MHz line. The only plausible molecule with 
a transition within 2 MHz of the observed freque ncy, accord- 
ing to the line catalogs of J. Cernicharo, CDMS (iMiiller et all 
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I2005h . and JPL (iPickett et all Il998l) . is '^CCH, whose Nj^f.^f 
- 47/2,4,7/2-37/2,4,7/2 transition lies at 336775.7 MHz. This 
molecule, however, is ruled out since the nearby 47/2,4,9/2- 
37/2,4,9/2 transition at 336756.2 MHz, with a slightly higher in- 
trinsic strength, is not present in our data (see Fig. [T]i. Since no 
other plausible candidate can be found for the 224525.4 MHz 
line and since unidentified lines of that intensity are rare in 
IRC +10216 at these frequencies, we conclude that we have al- 
most certainly detected CN . Confirmation of this identification 
would be highly desirable, but may not be easy to obtain. The 
next two rotational transitions of CN , at 449 and 561 GHz, can- 
not be observed from ground owing to high atmospheric opacity, 
and still higher J transitions may be too weak to detect in a cool 
source such as the outer envelope of IRC -1-10216. 

The J - 3-2 line of CN , which appears free of contam- 
ination by background lines, has a pronounced U-shaped pro- 
file, which for a spherical expanding envelope indicates that 
the emission is more extended than the half-power beam of the 
telescope (7" at 336 GHz). Thus CN" appears confined to the 
same outer envelope of IRC h- 102 16, as that of other molecu- 
lar anions observed in this source (e.g. [C ernicharo et al. 2007; 
iThaddeus et al.ll20 08: Cernicharo et al. 2008). A column density 
of 5 X 10'^ cm"^ and a rotation temperature of 16 K were de- 
rived from a rotational diagram constructed with the velocity in- 
tegrated intensities of the J = 2-1 and 3-2 lines given in Table[T] 
on the assumption of a uniform source with a radius of 20" which 
is typical of molecules distributed in the outer shell. The rotation 
temperature is consistent with CN" emission from the cool outer 
envelope. With a column density of the CN radical of 2 x 10'^ 
cm"^, derived from several hyperfine components of the A^ = 1-0 
and A^ = 3-2 transitions, we estimate a CN"/CN abundance ratio 
of 0.25 %, whi ch is comparable to th e C3N/C3N ratio in this 
source (0.52 %: IThaddeus et al.ll2008h . 

From the upper limit of the 7 = 2-1 line of C2H", a 3cr col- 
umn density of < 7 x 10'° cm"^ was derived on the assumption 
of a source with a radius of 20" and a rotation temperature of 
20 K. The estimated C2H"/C2H abundance ratio (< 0.0014 %) 
is at least 5 times smaller than the already small C4H/C4H ratio 
(lApundezll2009l) . 

3. Modeling and discussion 

To obtain a more reliable estimate of the abundance and excita- 
tion conditions of CN" in IRC + 1 02 1 6 we have carried out radia- 
tive transfer calculations based on the LVG formalism. The phys- 
ical pa rameters of the envelope have been taken from lA^unded 
(l2009l) . We included the first 20 rotational levels of CN". The 
rate coefficients for de-excitation by collisions with para-H2 
have been explicitly computed through quantum mechanical cal- 
culations for temperatures up to 70 K and for transitions involv- 
ing the first 9 rotational levels of CN". The calculations are de- 
scribed in the AppendixlAl For collisions with He, the rate coef- 
ficients computed for para-H2 were scaled down by a factor of 
1.37 (the ratio of the square roots of the reduced mass of each 
couple of collision partners). For transitions involving rotational 
levels higher than 7 = 8 the Infinite Order Sudden approxima- 
tion was used. As commented before, CN" is confined to the 
outer envelope of IRC +10216. We find that to reproduce the 
line profiles and relative intensities observed, the abundance of 
CN" relative to H2 must peak at a radius between 13" and 17" 
from the star. The adopted radial distribution, with a maximum 
abundance relative to H2 of 2.5 x 10"'^ reached at a radius of 15" 
(12" if expressed as particle density, see grey thick line in Fig.|2]), 
produces line profiles in reasonable agreement with the observed 
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Fig. 2. Abundance distribution derived for CN" in the envelope 
of IRC +10216 (thick grey Hne labeled as "CN" fit"), as it repro- 
duces the CN" observed line profiles (see Fig. [3]). Also shown 
are the abundances of CN", CN, and other molecular anions as 
calculated with the chemical model (multiplied by 0.0003, 5, 15, 
0.03, and 0.05 for CN, C2H", C4H", CeH", and C5N", respec- 
tively). The abundances are expressed as number of molecules 
per cubic centimeter. The angular distance is given in the top 
axis for an assumed distance to IRC +10216 of 120 pc. 




112260 112270 224510 224530 336760 336800 
Rest frequency (MHz) 

Fig. 3. Line profiles calculated with the LVG model (thick grey 
lines) using the compact CN" abundance distribution (thick grey 
line in Fig|2l) are compared with the observed CN" lines (black 
histograms). Fits to the CeH and SiC2 v^-l lines have been sub- 
tracted in the 7=1-0 and 2-1 observed spectra. The 7=1-0 line 
profile is very uncertain due to the blend with the strong CeH 
line. 



ones (see Fig.[3]l. We note that since the density decreases as the 
radius increases, the maximum in the particle density is reached 
at smaller radii than the maximum in the abundance relative to 
H2. The total column density across the envelope (twice the ra- 
dial value) is 3 x 10'^ cm"^, in good agreement with the value 
derived from the rotational diagram. In the region where most of 
CN" is present (at a radius of ~ 2 x lO'^ cm, where the gas ki- 
netic temperature is ~ 40 K and the density of H2 molecules 
is around 4 x lO'* cm"-*) the rotational levels involved in the 
CN" observed transitions are subthermally excited. Therefore, 
the collision rate coefficients utilized turn out to be essential to 
correctly estimate the CN" abundance in the outer layers of IRC 
+ 1 02 1 6 ' s envelope . 

To gain some insight into the formation of CN" in the ex- 
ternal layers of the molecular envelope of IRC +10216, we 
have perform ed chemical mode li ng cal culations similar to those 
described by ICernicharo et al.l (l2008h . The physical parame - 



ters of the envelope have been taken from ' Agunded (l2009l) . 
The rate constants and branching ratios of the reactions of 
ani ons with H, O, and N atoms, studied in the laboratory 
by lEichelberger et alJ (l2007h . have been updated according to 
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the values used by ICordiner & Milld (l2009t) and IWalsh et al] 
(20091J. Photodetachment rates of molecular anions have been 
assumed by Millar et al. ( 2007) to depend on the electron affinity 
of the neutral counterpart. For CN we have assumed the same 
rate expression adopted for CeH , as the neutral counterparts of 
both molecules have similar electron affinities (3.862 and 3.809 
eV respectivelv:'Rienstra-Kirac ofe et al.ll2002h . Plotted in Fig.|2] 
is the calculated radial distribution of the abundance of CN" 
(black thin line) and of some other molecular anions. CN" is pre- 
dicted to form at a much greater radius than C4H", CeH", C3N", 
and C5N", because, unlike the other anions, it is not formed di- 
rectly from the radical CN but through the reacti ons of the anions 
Q- (n = 5-10) with N atoms (see also lCordine'r & Millar 2009 b. 
CN being a small molecule, the rate constant for the reaction of 
radiative electron attachment is likely to be very small. Here we 
have assumed a value of 2 X 10~'^ cm^ s~^ at 300 K , similar to 
that computed for C2H by Herb st & Osamura (l2008l) . This pro- 
cess results in a too low formation rate for CN~, more than 5 
orders of magnitude lower than that provided by the reactions of 
C^ and N atoms. The reaction of HCN and H" is also a source 
of CN" in the inner regions of the envelope, but has only a mi- 
nor contribution (less than 0.2 %) to the total amount of CN" 
formed in the envelope. The anion C2H", on the other hand, is 
solely formed through the reaction of C2H2 and H", which takes 
place in the inner regions. According to our chemical model, 
CN" reaches a maximum abundance relative to H2 of 1.6 x 10"** 
at a radius of 8 x 10'^ cm, and a total column density across the 
envelope of 8 x 10'^ cm"^. For C2H", the model predicts a fairly 
small column density of 7 x 10'° cm"^, distributed over a region 
inner to 10"' cm. These results are in agreement with the recent 
chemical model of Cordiner & Millar ( 200 9), who predicted that 
both CN" and C2H" could be detected in the circumstellar enve- 
lope of IRC -Hi 02 16. 

The abundance and column density predicted for CN" by 
the chemical model is in reasonable agreement with the value 
derived from the observed lines and the LVG model. However, 
the calculated spatial distribution is markedly different from that 
derived by the observations (see Fig.|2l). In fact, by adopting the 
CN" abundance distribution obtained with the chemical model 
the resulting line profiles show important discrepancies as com- 
pared with the observed ones. While the calculated absolute 
line intensities are about the same order of magnitude as those 
observed, significant disagreements between the relative inten- 
sities and the line profiles are found. The calculated line in- 
tensity decreases too rapidly when going from the / = 1-0 to 
the J - 3-2, and the computed line profiles are too much U- 
shaped, with nearly all the emission predicted to occur at the 
line edges (i.e. at the terminal expansion velocity). These dis- 
crepancies arise because the chemical model predicts that CN" 
is present in a region of the circumstellar envelope that is too 
far from the central star. An abundance distribution more com- 
pact than predicted by our chemical model may arise if the en- 
velope is not modeled as being smooth, but as having density- 
enhanced shells. Cordiner & Millar (2009) have recently studied 
the effect of such density enhancements on the radial distribu- 
tion of molecular abundances and found that molecules formed 
in the outer envelope would concentrate at the position of the 
first and/or second shells, located at 15 and 27" respectively. 

For C2H", the chemical model predicts it to be distributed 
over an 8" diameter region (see Fig.O with a total column den- 
sity of 7 X 10'° cm"^. Once averaged over the 14.6" beam of the 
IRAM 30-m telescope at the frequency of the / = 2-1 transition. 
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the calculated column density is about 3 times lower than, and 
thus consistent with, the 3cr upper limit derived from the non 
detection of the 7 = 2-1 line. 

The identification of CN" in IRC -1-10216 with a relatively 
large anion-to-neutral abundance ratio (0.25 %) suggests that it 
may be detectable in other astronomical sources. Upper limits to 
the CN"/CN abundance ratio as low as 0.2-2 % were obtained 
in TMC-1, L1527, Barnard 1, and the Orion Bar in a previous 
search for the J = 2-1 transition by Agiindez et al. (2008). More 
sensitive observations would be needed if the abundance of CN" 
in other sources is similar to that found in IRC +10216. 

The high abundance of CN" compared to that of C2H" 
demonstrates the efficiency of the reactions of N atoms and large 
carbon anions. A more sensitive search for C2H" might support 
this alternate scheme for the formation of anions in space, and 
perhaps explain the low observed abundance of C4H" as com- 
pared to C3N". 
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Appendix A: CN -H2 collision rate coefficients 

The potential energy surface (PES) of the CN -H2 complex was 
calculated ab initio using single and double-excitation coupled 
cluster method with non-iterative triple excitations [CCSD(T)] 
( JKnowles et al.Lll993l l2000) implemented in MOLPRcB The ge- 
ometry of the system was described in the body-fixed frame 
and characterized by three angles {6, 0', <p) and the distance 
R between the centers of mass of H2 and CN . The H2 bond 
distance was fixed at ro= 1.44876 ao and the CN" bond dis- 
tance was varied for the purpose of the averaging of the PES 
over the lowest vibrational state of the CN" diatom. The ba- 
sis set superposit i on err or correction counterpoise procedure of 
iBovs & Bemardil (11970 ) was applied. The four atoms were de- 
scribed by th e correlati on-consis tent tr iple zeta basis set (aug- 
cc-pVTZ) of IWoon & Dunning (i 19941) augmented by t he (3s , 
2p, Id) midbond functions defined by IWilliams et al.l (1 19951) . 
placed at mid-distance between the CN" and H2 centers of mass. 
The final V(r,R,d,ff ,(p) PES is five-dimensional, however, in 
this work we included only three perpendicular orientations of 
the H2 molecule {{ff , 0) pairs: (0, 0), (0, 90), (90, 90)] to average 
over H2 rotations. Additionally, the PES was averaged over the 
CN" internuclear distance corresponding to the CN" vibrational 
ground state wave function. The 2-D PES was finally obtained 
as an arithmetic average of three H2 orientations. The full five- 
dimensional PES and four-dimensional scattering calculations 
will be presented elsewhere. 

We considered collisions of CN" with para-H2(7'2 - 0) at 
low temperatures. The rotational levels of CN" and H2 are des- 
ignated by j\ and J2, respec tively. We used the f u lly qu antum 
close-coupling approach of [Arthurs & Dalgarna (Il960l) . The 
standard time-independent coupled scattering equations were 
solved using the MOLSCAT code (Hutson & Green, 1994. 
Calculations were carried out at values of the total energy rang- 
ing from 3.6 to 500 cm"'. The integration parameters were cho- 
sen to ensure convergence of the cross sections over this range. 
At the largest total energy considered (500 cm"') the CN" ro- 
tational basis included channels up to y'l = 21 to ensure con- 
vergence of the excitation functions a-j^p^f j'^iEc) for transitions 
including up to the ji - 8 rotational level of CN" . The rotational 
basis of H2 was restricted to 72 = levels. The coupling with the 
72 - 2 (a nd higher) states of H2 was not taken into account. As 
shown bv lLique et al.l (l2008h . this approach is expected to yield 
reUable results for the energy range considered here. From the 
above described excitation functions one can obtain the corre- 
sponding state-resolved thermal rate coefficients by Boltzmann 
averaging: 
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Fig.A.l. Collisional de-excitation rate coefficients of CN" by 
para-H2 are shown as a function of temperature for the J - 1-0, 
2-1, 2-0, and 3-1 rotational transitions of CN". 



Table A.l. CN -H2 


collision rate coefficients (10 ' cm-' s ') 








Temperature 


(K) 






Transition 


10 


20 


30 


40 


50 


60 


70 


1 ^0 


3.33 


3.06 


2.88 


2.72 


2.59 


2.48 


2.39 


2^0 


0.81 


0.65 


0.59 


0.55 


0.54 


0.53 


0.53 


2^ 1 


4.81 


4.50 


4.25 


4.05 


3.88 


3.74 


3.62 


3^0 


0.67 


0.64 


0.61 


0.59 


0.57 


0.56 


0.52 


3^1 


1.32 


1.18 


1.07 


1.00 


0.96 


0.93 


0.92 


3^2 


4.81 


4.62 


4.39 


4.19 


4.03 


3.88 


3.77 



Note.- The complete set of de-excitation rate coefficients of CN" by 
collisions with para-H2 considered in this study is available on the 
BASECOL website http : //basecol . obspm . f r/ 



coefficients up to 70 K. The temperature dependence of the rate 
coefficients for selected de-excitation transitions is illustrated in 
Fig. lA.ll with the values given in Table I A. fl 
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where k is the Boltzmann constant. To obtain precise values of 
the rate constants, the energy grid was chosen to be sufficiently 
fine to include the numerous scattering resonances. The total en- 
ergy range considered in this work allows us to determine rate 



^ MOLPRO, version 2006.1, a package of ab initio programs, H.-J. 
Werner, P. J. Knowles, R. Lindh, F. R. Manby, M. Schutz, R Celani, T. 
Korona, G. Rauhut, R. D. Amos, A. Bernhardsson, A. Berning, D. L. 
Cooper, M. J. O. Deegan, A. J. Dobbyn, F. Eckert, C. Hampel and G. 
Hetzer, A. W. Lloyd, S. J. McNichoIas, W. Meyer and M. E. Mura, A. 
Nicklass, P. Palmieri, R. Pitzer, U. Schumann, H. Stoll, A. J. Stone, R. 
Tarroni and T. Thorsteinsson, see http : //www . molpro . net 



